1 



O 
O 

(N 

P-h 



> 

o 

(N 
O 
(N 
O 

">< 

Oh! 
<D . 
43 



X 



A Measurement of the Charge Asymmetry 

(Submitted to Physical Review Letters, December 31st, 2001) 

A. Alavi-Harati 12 , T. Alexopoulos 12 , M. Arenton 11 , K. Arisaka 2 , S. Averitte 10 , R.F. Barbosa 7 '**, A.R. Barker 5 
M. Barrio 4 , L. Bellantoni 7 , A. Bellavance 9 , J. Belz 10 , R. Ben-David 7 , D.R. Bergman 10 , E. Blucher 4 , G.J. Bock 7 
C. Bown 4 , S. Bright 4 , E. Cheu 1 , S. Childress 7 , R. Coleman 7 , M.D. Corcoran 9 , G. Corti 11 , B. Cox 11 , M.B. Crisler 7 
A.R. Erwin 12 , R. Ford 7 , A. Glazov 4 , A. Golossanov 11 , G. Graham 4 , J. Graham 4 , K. Hagan 11 , E. Halkiadakis 10 
J. Hamm 1 , K. Hanagaki 8 , S. Hidaka 8 , Y.B. Hsiung 7 , V. Jejer 11 , D.A. Jensen 7 , R. Kessler 4 , H.G.E. Kobrak 3 , J. LaDue 5 
A. Lath 10 , A. Ledovskoy 11 , P.L. McBride 7 , P. Mikelsons 5 , E. Monnier 4 **, T. Nakaya 7 , K.S. Nelson 11 , H. Nguyen 7 ** 
V. O'Dell 7 , M. Pang 7 , R. Pordes 7 , V. Prasad 4 , X.R. Qi 7 , B. Quinn 4 ^, E.J. Ramberg 7 , R.E. Ray 7 , A. Roodman 4 
M. Sadamoto 8 , S. Schnetzer 10 , K. Senyo 8 , P. Shanahan 7 , P.S. Shawhan 4 , J. Shields 11 , W. Slater 2 , N. Solomey 4 
S.V. Somalwar 10 , R.L. Stone 10 , E.C. Swallow 4 ' 6 , S.A. Taegar 1 , R.J. Tesarek 10 , G.B. Thomson 10 , P.A. Toale 5 
A. Tripathi 2 , R. Tschirhart 7 , S.E. Turner 2 , Y.W. Wah 4 , J. Wang 1 , H.B. White 7 , J. Whitmore 7 , B. Winstein 4 
R. Winston 4 , T. Yamanaka 8 , E.D. Zimmerman 4 

The KTeV Collaboration 

1 University of Arizona, Tucson, Arizona 85721 
2 University of California at Los Angeles, Los Angeles, California 90095 
3 University of California at San Diego, La Jolla, California 92093 
4 The Enrico Fermi Institute, The University of Chicago, Chicago, Illinois 60637 
5 University of Colorado, Boulder, Colorado 80309 
6 Elmhurst College, Elmhurst, Illinois 60126 
7 Fermi National Accelerator Laboratory, Batavia, Illinois 60510 
8 Osaka University, Toyonaka, Osaka 560-0043 Japan 

9 Rice University, Houston, Texas 77005 
10 Rutgers University, Piscataway, New Jersey 08854 
11 The Department of Physics and Institute of Nuclear and Particle Physics, 
University of Virginia, Charlottesville, Virginia 22901 
12 University of Wisconsin, Madison, Wisconsin 53706 



We present a measurement of the charge asymmetry Sl in the mode Kj, 



z e T v based on 298 million analyzed 



decays. We measure a value of Sl = (3322 ± 58(stat) ± A7(sys)) ■ 10~ 6 , in good agreement with previous 
measurements and 2.4 times more precise than the current best published result. The result is used to place more 
stringent limits on CPT and AS = AQ violation in the neutral kaon system. 

PACS numbers: ll.3Q.Er, ll.30.-j, 13.25.Es, 13.20.-v 



The charge asymmetry in Kl semileptonic decays is 
deeply related to CP violation through neutral kaon 
mixing. The Kl wavefunction is proportional to (1 + 
el)K — (1 — eL)K°, where €l parameterizes CP vio- 
lation in mixing. This parameter can be determined in 
the decay modes Kl — > 7r ± e =F ^ (Ke3) and Kl — > 7r ± /i T ^ 
(K/x3). The charge asymmetry (Sl) is defined as 



Sl = 



BR(e 



- BR(e~7r+) 



BR(e-ir+) 



Assuming AS = AQ and no CPT violation in the decay 
K°(K°) -> e + n-v(e'n + v), S L is simply 2Rc e L . To first 
order in the parameters that violate CP and CPT, 

S L = 2Re e L - 2Re Y - Re(x - x). 

The terms Y and Re(x — x) parameterize CPT viola- 
tion in the AS — AQ and AS = —AQ transitions 
respectively]^, S: 

< e+TT-v\K° > l-Y 



< e-n+D\K a >* 



1 + Y 



< e+7r v\K° > 



< e-TT+9\K° >* 



In the Standard Model, the AS" = AQ violation is CPT- 
conserving and occurs in 2nd order weak decays. Esti- 
mates for |jc| are in the range 10 _7 ||. Sl is modified to 



order 



and will not be considered in this letter. 



The measurement of Sl can be compared to expec- 
tations derived from the Kl — » 2-7T (K7r2) amplitudes 
(77+_, 7700) 1,1: 



Re (Y 



+ a) 
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This comparison is sensitive to CPT violation in Ke3 and 
K7r2 decays, where the latter is parameterized by Re a. 

The current world average on Sl comes mainly from 
a measurement by the CERN-Heidelberg collaboration 
in 1974. Using 34 million Ke3 decays, they found 
S L = (3409 ± 171(stat) ± 50(sys)) ■ 10- 6 @. We 
describe in this letter a measurement of Sl in the Ke3 
mode based on approximately 298 million Ke3 decays. 
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FIG. 1: The KTeV detector configured for measuring 
Re(e'/e). 



The KTeV beamline and detector at Fermilab is de- 
scribed elsewhere ||. For this analysis, the detector is 
configured for the measurement of Re(t'/e). As shown 
in Fig. [lj two approximately parallel neutral Kl beams 
enter a long vacuum tank, which defines the fiducial vol- 
ume for accepted decays. One of the beams strikes an 
active absorber (regenerator), which serves to tag the co- 
herent regeneration of Kg- The regenerator moves to 
the other beam in between Tevatron spill cycles. Follow- 
ing the vacuum tank, there are 4 planar drift chambers 
and an analysis magnet that imparts a 411 MeV/c hor- 
izontal transverse kick to the charged particles. A high 
precision 3100-element pure Cesium Iodide calorimeter 
(Csl) is used primarily to measure the energy of e and 
photons. Photon veto detectors surrounding the vacuum 
tank, drift chambers, and Csl serve to reject events with 
particles escaping the calorimeter. 

Only Ke3 decays in the beam opposite the regenera- 
tor were used in this analysis. They were triggered by 
the presence of hits in the drift chambers and scintil- 
lators placed immediately upstream of the Csl, and by 
the lack of activity in the regenerator or photon vetoes. 
To keep the event rate managable, the trigger rejected 
beam muons and K/z3 events by requiring no activity in 
the scintillation counters (muon veto) placed behind 4 
meters of Fe absorbers downstream of the Csl. 

Selecting Ke3 events is relatively straightforward, as 
they are the most copious Kl decay mode (BR(Ke3) 
« 0.39). We identify two-track vertices in the region 
between 90 and 160 meters from the target with kine- 
matics inconsistent with Kl — ► 7r + 7r~7r° (K7r3) and 
A(A) — > p7r - (p7r + ). The tracks are required to extrap- 
olate to Csl energy deposits while maintaining sufficient 
clearance from the inner and outer edges. A track is iden- 
tified as an electron if its momentum exceeds 5 GeV/c 
and its Csl energy deposit divided by its momentum 
(E/P) exceeded 0.925. Pions are identified as having mo- 
mentum exceeding 8 GeV/c and E/P less than 0.925. 
The tighter pion momentum requirement rejects K^3 



background, since 8 GeV/c is the threshold for minimum- 
ionizing particles to penetrate through to the muon veto 
counters. 

To minimize background from target Ks coherently 
interfering with the Kl , the event proper time (r) had to 
exceed 10.5 Ks lifetimes (rg). The Kl momentum (Pk) 
reconstruction has ambiguities inherent to any 3-body 
decay with one unobserved particle. The longitudinal 
component of the neutrino momentum in the kaon rest 
frame can be either parallel or anti-parallel to the kaon 
flight direction in the lab frame. This introduces two 
possible solutions for Pk- We use the low Pk solution 
to calculate r. Based on Monte Carlo (MC) simulations, 
this choice is at least 70% correct for r < 22.5 rg. A small 
correction, derived from MC, was made for errors in the 
Pk solution and the residual target Ks interference. 

The measurement of 8l requires a careful control of 
systematics. The most important of these is to ensure 
an equal acceptance and efficiency for oppositely charged 
particles. We combine data sets of opposite analysis mag- 
net polarities, which was reversed about once per day, to 
ensure that the detector is exposed to oppositely charged 
particles in the same manner. We also used data collected 
during rare decay running (E799), which included a tran- 
sition radiation detector (TRD)|?], |) used for the study 
of electron and pion identification. 

There are eight possible configurations for Ke3 decays: 
e + 7r~ or e~7r + , east (E) or west (W) Kl beams, and pos- 
itive (+) or negative (— ) magnet polarities. The number 
of Ke3 events in each configuration (iVj) depends on the 
branching ratio (BR), the acceptance and efficiency (A), 
and the flux N(Kl)- We define R as: 

4 = BR(e+n-)A(e+7r-, E, +)N(K L ,E, +) 

BR{e-<K+)A(e-%+,E,-)N(K L ,E,-) [ ' 
BR{e+TT-)A(e+Tr~,E, -)N(K L , E, -) 
BR{e-Tr+)A(e~Tr+, E, +)N{K L , E, +) 
BR(e+7r-)A(e + Tr-, W, +)N(K L , W, +) 
X BR(e-7r+)A(e~Tr+, W, -)N(K L , W, -) 
BR(e+TT-)A(e+Tr-,W, -)N(K L , W, -) 
BR{e-Tr+)A(e~Tr+, W, +)N(K L , W, +) ' 

where the four numerators and denominators represent 
Ni, i = 1,8. Since the fluxes cancel and Afe^Tf^) = 
A(e T ir ± ) under magnetic field reversal, Eq. || reduces 
to: 



raw Sl = 



R-l 



1 



a(mw8 L ) = --,mtiTF- ( 3 ) 



1 



R+l v ' 8 V Ni 

Equation ^ defines the "raw" Sl, and does not take into 
account several effects discussed below. Our data yields 
raw(5 L = (3417 ± 58) • 10~ 6 (ppm). Figure | shows its 
dependence on the decay vertex distance from the target. 

In principle, the same flux and acceptance/efficiency 
cancellation occurs in the pairing of the first and second 
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FIG. 2: Raw Sl and number of events versus vertex distance 
from target. A fit to a constant yields 3417 ± 58 ppm with a 
X 2 of 81.53 for 81 degrees of freedom. 

(third and fourth) factors of Eq. 0, such that Sl can be 
measured within the east and west beams. Eq. |^ is the 
geometric mean of the two Sl measurements. 

Figure shows the raw asymmetry in various data 
subsets. In particular, we see the importance of com- 
bining data from opposite polarity runs. The difference 
between measurements made with only one magnet set- 
ting, <5l(+) — Sl(-), is significant and on the order of 
Sl itself. This difference is due to very small offsets of 
the inner apertures and is well- reproduced by our MC. 
The MC is not used to correct this large effect since 
Eq. H naturally accounts for such geometrical asymme- 
tries. On the other hand, the values for Sl agree between 
decays in the east and west beams. Other data divi- 
sions yield good agreement except for the sample with 
8 GeV/c < < 15 GeV/c (discussed below). 
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FIG. 3: Raw Sl for all events and various data subsets. The 
line indicates the raw value of 3417 ppm. 

We account for the different behavior of particles and 



antiparticles in matter. For each cut that introduces a 
bias, we measure / + (f~), the inefficiency in the e + ir~ 
(e~7r + ) configuration. The correction to raw Sl for the 
cut, summarized in table |, is simply (/ + — f~)/2. 

Many biases due to ir^ interaction differences were con- 
sidered, and data were used to measure most of these 
corrections. These include the energy deposition in 
the Csl, the loss of pions due to interactions in the trig- 
ger scintillators, and pion punchthrough past the Fe ab- 
sorbers depositing energy in the muon veto. Events with 
E/P for pions exceeding 0.925 would be removed since the 
analysis requires exactly one identified electron. This ef- 
fect, studied using pions in K7r3 events, removes ~ 0.6% 
of the pions, with the probability being slightly larger for 
7r + than 7r~. The asymmetry between 7r + and 7r~ has 
a strong momentum dependence and explains the trend 
seen for the 8 GeV/c < P v < 15 GeV/c data sample (see 
Fig. ||). The correction to Sl is — 156 ± 10 ppm. 

A correction was estimated for possible biases due 
to the trigger muon veto requirement, which removed 
approximately 3% of pions due to decay-in-flight and 
punchthrough. Good events can also be lost if accom- 
panied by accidental muons. Only pion punchthrough 
would cause a bias due to the neutron excess in the Fe 
absorbers. The estimate used a Ke3 sample collected in 
E799, where the TRD gave additional electron purity. 
The sample is triggered by requiring muon counter hits, 
as these events are the ones that would be removed by the 
analysis trigger. This sample had a charge asymmetry of 
(4.4 ± 1.3) • 10 -3 , consistent with no bias. Nevertheless, 
we use a correction of 34 ± 40 ppm, where the error is 
given by the sample statistics. 

Pion interactions in the trigger scintillator can confuse 
the reconstruction so that the track fails to match to a 
Csl cluster. This effect, studied using Ke3 events with 
a relaxed track-cluster matching requirement, removes « 
0.42% of all pions and was slightly more probable for the 
7r~ than 7r + , leading to a correction of 54 ± 10 ppm. 

Considering interaction differences, one has primar- 
ily e + annihilation and <5-ray production. An e + annihi- 
lation occuring near the upstream surface of the trigger 
scintillator would fail the trigger requirement due to the 
reduced energy deposition. <5-rays more than 10 MeV and 
emitted from the vacuum window or tracking system can 
cause losses since the momentum transfer would drasti- 
cally change the event kinematics. These losses cause 
a bias due to the small difference between the Bhabha 
(e + e~) and Moller (e~e~) scattering cross section. These 
corrections were derived by a Geant simulation of the de- 
tector material. They are of order ±10 ppm, mitigated 
by the high momentum of the tracks and the minimal 
material (2% Xq) upstream of the Csl. 

Finally, a Ke3 sample was used to study potential bi- 
ases of the E/P cut on electrons. This sample was se- 
lected by a relaxed E/P cut and by the tagging of a 
minimum- ionizing 7r in the Csl. The E/P removed re- 
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moved approximately 0.233% of e + and e~. No bias is 
seen, and we assigned a correction of —19 ±18 ppm. This 
result was confirmed with a similar Ke3 sample collected 
in E799, where the TRD gave additional electron purity. 

We applied a small correction to account for pions ab- 
sorbed in the vacuum window, drift chambers, and he- 
lium bags. This was estimated assuming isospin con- 
servation in strong interactions, and so the 7r + and tt~ 
absorption differences depended on the excess protons, 
which was predominantly in the form of hydrogen. We 
accounted for a small bias due to ir~ charge exchange 
near the upstream surface of the trigger scintillator (anal- 
ogous to the case of e + trigger loss). For this, we used 
the measurements of || and JjiJ and extrapolated them 
to our 7r~ momentum range. 

Other small corrections, derived from MC, accounted 
for the inexact reversal of the analysis magnet polarity 
and the residual coherent Ks from the target, absorber, 
and Ki, scatters with the final collimator and regenera- 
tor. Table | shows the summary of all systematic cor- 

TABLE I: Corrections in ppm for this analysis. 



tv + it~ difference in Csl 


-156 ± 10 


tt ± interaction in trigger scintillator 


+54 ± 10 


7r decay and punchthrough 


+34 ± 40 


e + e~ difference in Csl 


-19 + 18 


target/absorber Ks — Kl interference 


-12 ± 1 


e + annihilation in spectrometer 


+ 11 ± 1 


Delta ray production 


-8.5 + 4.3 


7r absorption in spectrometer 


+5.0 + 3.2 


inexact analysis magnet polarity reversal 


-3.1 ± 1.6 


final collimator and regenerator scatters 


-1.2 + 2.3 


■K + ir~TT° background 


+0.5 + 0.7 


K/j.3 background 


+ 


A p7T and background 


+ 


Total correction 


-95.3 + 46.5 



rections. The uncertainties are uncorrelated since they 
are mainly statistical errors of the Ke3 and K7r3 control 
samples. Combining the corrections with raw 5l, we find: 

5 L = 3322 ± 58 (stat) ± 47 (syst) ppm 
= 3322 ± 74 (combined) ppm. 

This result is in excellent agreement with previous mea- 
surements and 2.4 times more precise than the current 
best result (see Fig. |]). A combination of all results 
including ours yields: 



5 l = 3307 ± 63 ppm 



4.2/6 d.o.f. 



Substituting 77^ and 7700 values from and the com- 
bined 5l into Eq. [I], we find: 



The result is consistent with no CPT violation and is 
limited by the charge asymmetry uncertainty. It limits 
\Re (Y + (x - x)/2 + a)\ < 61 ppm at the 90% C.L.. 
Barring fortuitous cancellations, this is thus far the most 
stringent limit on Y, x, and a. 
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FIG. 4: Compilation of Sl measurements including this result. 
Lines indicate the average and its uncertainty. 
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Re (Y 



a) = (1650+ 16) - (1653 + 32) 
= —3 ± 35 ppm. 
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